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g . 1 Introduction 



This work was motivated by the paper of Shkalikov |4]] concerning the analysis of 
the semi-classical Airy operator, together with our computer simulations of the 
associated discrete problem using the numerical package Matlab. Specifically, we 
examined the operator H$ h given formally by 

H s , h :=-h 2 -^ + V 5 onL 2 (-l,l) (1) 

where 

J i(x + 5) for x > 



V s (x) 



i(x — 5) for x < 



and both 5 > 0, h > are small. In it is shown rigorously that as h — > the 
spectrum of Hq ^ becomes dense inside an arbitrarily small neighbourhood of the 
Y-shaped subset of C denned by 

[z,l/V3],[-z,l/V3] and [1/^3, oo), 

where we use [c£,/3] to denote the line segment joining a, (3 G C (see figure 2). 
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This paper confirms our surprising numerical results, which suggested that when 
an arbitrarily small jump discontinuity is inserted into the otherwise linear (purely 
imaginary) potential, and then the semi-classical parameter h allowed to go to 0, 
the asymptotic spectrum of Hs,h turns out to be completely different from that of 
Ho,h (see figure 1 and Corollary |). 

Several papers [§, ||, U have been written about the operator H h - a major moti- 
vation being that it is a model operator for the 'Orr-Sommerfeld' problem The 
operator also defines the 'Squire model for the Couette flow' in hydrodynamics; 
and in its own right, defines the semigroup which is the solution of the so-called 
'Torrey equation' ||, related to the diffusion of magnetic fields. Thus, although the 
spectrum of this non- self- adjoint operator displays sometimes strange and singular 
behaviour, it must not be dismissed as a 'pathological' example from pure math- 
ematics since it has important applications - for example, in magnetic resonance 
imaging devices. 

It is well known that a basis for solutions of the so-called 'Airy equation' 

-f"(z) + zf(z) = 

is given by any two of the Airy functions Ai(z), Ai(e~ 2n ^ 3 z) and Ai(e 2n ^ 3 z) (see 
||). Thus, the analytic investigation of the eigenvalues will involve examining the 
asymptotic behaviour of certain Airy functions, where the expressions we use are 
WKB-type approximations. We will show that the eigenvalues lie inside a certain 
subset of the complex plane, which is intimately related to the Stokes' lines (or 
principal curves) of the problem (see ||). The proof will depend upon showing 
that for all A outside this subset, the eigenvalue problem 

H h f(x) = Xf(x) 

has a well-defined Green's function. As in [|J , our analysis uses the concept of the 
characteristic determinant which we will describe in Section 

In Section |4] we analyse the spectral behaviour in the semi-classical limit h — > for 
a general complex-valued, piecewise linear potential. The surprising result is that 
each linear segment of the potential gives rise to a characteristic 'Y-shaped' set of 
eigenvalues; and the spectrum of the operator is contained within the superposition 
of these 'Y-shaped' sets. From the point of view of applications, this means that 
whilst the asymptotic spectrum is theoretically computable for an idealised linear 
potential; in practice, any arbitrarily small deviation from the ideal can completely 
change the spectrum. One way of expressing this for an operator H is in terms of 
the pseudospectral sets: 

Spec e (#) := Spec(#) U {z G C : \\{H - z)- l \\ > e" 1 }, 

i.e. the contour sets of the resolvent norm, with the convention that z G Spec(if ) 
implies 

\\(H - z)- l \\ := oo. 
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For many non-self-adjoint operators it has been demonstrated |3|, |^| that the 
pseudospectral sets become very large as some parameter varies, even though z 
may be far from the spectrum of the operator. This is equivalent to saying that the 
spectrum is computationally unstable. Our aim in this paper is to demonstrate for 
a relatively transparent case, the mechanism behind this phenomenon; we believe 
the results to be capable of extension to a more general class of piecewise analytic 
potentials. 

In Section |^ we provide an analysis of the simultaneous limit as h — > and 5 — > 
together. 



2 The Characteristic Determinant 



In this section we describe the characteristic determinant of the operator defined 
by 



H h f(x) := -K- 



dx 2 



+ V{x)f{x) 



acting on L 2 (— 1, 1) with Dirichlet boundary conditions, h > small, and V{x) the 
complex valued, ra-times piecewise linear function 



V(x) 



m\x + l\ xo < x < xi 

77122 +I2 X\ < X < X2 



k 777 n X ~\~ l n X n —\ <C X ^ X n 

with — 1 = x < x\ < . . . < x n = 1, and the m^, ij i = 1, . . . , n complex constants. 
The domain of the operator is given precisely by 

Dom(# fc ) = {fE C[-l, 1] : /(-l) = /(l) = 0, /' G C[-l, 1] and f" G L 2 (-l, 1)} 

(2) 

where the primes ' denote differentiation with respect to x, and /" is initially to 
be interpreted in the distributional sense. A direct substitution shows that a basis 
of solutions for the differential equation 

-h 2 f"(x) + (V(x) - X)f(x) = 

where V(x) := mx + I; and /, m are complex constants, is given by any two of the 
Airy functions Ai(w) and Ai(e ±2m ^ 3 w), where 

w ■= h- 2 ^mr 2 ' z {V{x) - A). 

It follows that, in order to construct an eigenfunction of the operator Hh, we seek 
constants an, i = 1, . . . , n such that the function 



/(*) == 



CtnUii(x) + Ol\2Ux2{x) X Q < X < Xi 

a2iu 2 i(x) + a 2 2U 2 2(x) x\ < x < x 2 
, a n iu n i(x) + a n2 u n2 (x) x n -i < x < x t 



satisfies all of the domain conditions (0), where 



u a (x) := Az(e- 2fcm/3 /i- 2/3 m 2 72/3 (K^ + h) ~ A)) (3) 
with k e { — 1,0,1}. For each i = 1, . . . , n, the functions Ui 2 are defined similarly, 
except that a different choice of k is to be taken from { — 1, 0, 1}. 

In addition to satisfying the boundary conditions /(— 1) = /(l) = 0, / must also be 
continuously differentiable, even at the points x«. From the power series definition 
0, p. 54], it is clear that the Airy functions Ai are analytic on the whole of C, and 
so the requirement that / be continuously differentiable reduces to the 2(n — 1) 
simultaneous 'matching' conditions 



atiiUix(xi-) + a i2 u i2 (xi-) - a>(i+i)iU( i+1 ) 1 (x i +) - a>(i+ipU(i + ip(xi+) = 



and 



anu'^Xi-) + a i2 u' i2 (xi-) - a { i +l)1 u[ i+t)l (x i +) - a { i +1)2 u[ i+l)2 (x i +) = 0. 
The boundary conditions /(— 1) = /(l) = demand that 

aii«ii(-l) + a 12 u 12 (-l) = 

and 

a n iM„i(l) + a n2 u n2 (l) = 0. 
Thus finding a solution of the differential equation 

-h 2 f"(x) + V(x)f(x) = Xf(x) 

which satisfies all of the domain conditions (0), involves solving the matrix equation 



(-1) 


u 12 (-l) 
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for the constants aji, a^- Note that we are taking the left- and right-hand limits at 
the points Xj, although here and subsequently we will abuse the notation in order 
to add clarity, and simply write Un(xi) etc. It is the determinant of the matrix 
on the left-hand side of @) that we shall call the characteristic determinant of the 
eigenvalue problem defined by Hh- 
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3 Airy functions and Stokes' lines 



For our proofs in the next section, we will need some notation and basic properties 
of the Airy functions. In all that follows we let the argument function Arg take 
principal values i.e. 

Arg : C — > (— 7r, 71"]. 

If Arg (j3 — a) := 9, the subsets Y(a, (3) of C are to be constructed as follows: take 
the lines 

a + re ™/a and f/3 + re 2 ^ 2 -/ 3 for 9 < 
a + TG and \p + re 26 */ 3 - 2 ™/ 3 for 9 > 

as r ranges in [0, oo), to their point of intersection, T say. Then, from T extend 
the infinite line defined by the set of A G C such that 

Re {(e- 29 ^ 3 (a - A)) 3/2 } = Re {(e- 29 ^((3 - A)) 3/2 }. (5) 

The motivation for this set will become clear during our proofs; in fact, it will be 
seen to comprise a curve asymptotic to the line 

L eC :I mM = ^M±Wg)}. 

Note for now, however, that (|5|) is h- independent. The ^-neighbourhood of any 
subset T of C will be defined by 

Nhd(T; e) := {t + z : t G T and \z\ < e}. 

We will use the well-known [Q] asymptotic expansion of the Airy function Ai, giving 
the WKB-type approximation: 

Ai(z) = exp (-\z 3 ' 2 ) (1 + 0{z-^)) (6) 

as \z\ — > oo, valid for |Arg (z)\ < 7r; and where the principal value of z 3 ^ 2 is taken. 
Following the notation of Olver (see p.413]), we define 

S := {z:\Arg (z)\<n/3} 

S 1 := {z : tt/3 < Arg (z) < n} 

£>_! := {z : -tt/3 > Arg (z) > -n} 

(suffixes enumerated modulo 3). One can check that for all complex z (and taking 
principal values), we have 



-Re {{e 2m l 3 zf' 2 } for z G & U 5_i 
Re {(e 2 ™' 3 z) 3 / 2 } for z G S 



Re { (e--/ 3 z) 3 / 2 } = , ' " i * ? U ^ (7) 



Re {(e~ 2 ^zf 2 } = I - Re n { [ Z £x } T * ^ c U S ° («) 

lv ; J [ Re {(z) 7 } for 2 G 5_i w 

and 

Ke 116 2j 1 _ \ Re {{zfl 2 } for z G S x . [J) 
Then, putting 

Ai k (z) := Ai(e- 2kni/3 z) (10) 

the asymptotics ([]) show that as |z| — > 00, |Azfc(z)| decreases exponentially for 
z G Sfc, and increases exponentially for 2 G S^-i U Sfc+i. The boundaries of the 
sectors i.e. the rays te ±m ^ 3 and te m for t G [0, 00), are known as the Stokes' 
lines (or principal curves) of the problem [0, p. 503]. Indeed, for the Airy equation 

-f"(z) + zf(z) = 

the Stokes' lines are defined to be the values of z such that 

2 



Re jT \Tt dt = Re {g^ 3/2 } = 0, 



and denote the boundaries of the principal subdomains S\ etc., inside of which 
the asymptotic expression (§) is valid for each k. 

We will call the suffix k 'allowable' for any given z G C, if 

Arg (e~ 2km/3 z) <tt. (11) 



4 The Semi-Classical Limit 

The following is a generalisation of the argument used in |4j] for the potential 
V(x) = ix, and will form a lemma for our main theorem. 

Proposition 1 Let V be the complex valued linear potential given by 

V(x) = mx + I x G [—1, 1] 

where m and I are complex constants; un(x), Uuix) are as defined in (^), and 
a,b G [—1, 1], a < b. 

Let e > be given and A G C. If 

\$Nhd(Y(V{a),V(b));e), 

then 

u u (b)u 12 (a) = o(uu{a)u 12 {b)) (12) 

as h — > 0. 
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Proof A simple scaling and translation of the operator Hh allows us, without loss 
of generality, to assume that a := — 1, b := 1 and 1 = 0. That is, we assume 

V(x) := xe* e , where 9 := Arg (m). 

By elementary trigonometry, one can check that we then have 



4 (\e\ 

—= sin — 

v/3 U 



e" + -= sin I — I e 2(e ~ T)l/3 



or 

Recalling our definition of the Airy functions un(x) and tt^Oz) (H), we put 

z{h,X,x) :=h~ 2 / 3 e- 2et / 3 (xe ei -\), 
and can rewrite (f|) explicitly in terms of h. Then, taking the modulus we obtain 

\Ai(z(h, X,x))\ 

ei ^ - 1 / 4 

= /i 1/6 =- exp -4Re (e- 2 ^ 3 (xe^ - A)) 3/2 ] (1 + 0(h)) (13) 

2^71" V 3h ) 

as h — >• 0, valid for |Arg A, < 7r. Therefore, in order to estimate the 
moduli of the Airy functions Aik(z(h, A, x)) in the limit h — > 0, it is sufficient to 
examine the behaviour of the functions 

x i-h. Re {(e- 2fcm / 3 z(/i, A, x)) 3/2 }, x G R 

for A; = -1,0,1. 

The basic idea of our proof is to show that for A outside an arbitrarily small e- 
neighbourhood of Y(V(— 1), ^(1)), one can assign allowable values of j and k from 
{ — 1, 0, 1} (in the sense of ([TT|) ) such that the inequalities 

Re {(e- 2jni/3 z(h } A, -1) 3/2 } < Re {( e - 2jm/3 z(h, A, 1) 3/2 } (14) 

and 

Re {( e - 2kni / 3 z(h, A, 1) 3/2 } < Re {(e- 2h7Ti/3 z(h, A, -l) 3 / 2 } (15) 

hold in the limit h — » 0. This will then be enough, by ([13]), to ensure that 
uu(l)u 12 (— 1) and nn(— 1)mi 2 (1) are of different orders of magnitude as h — » 0, 
thus implying (|l"2j). 

Using the statements of the previous section; for all values of A such that 

z(h,X,±l) :=h- 2 / 3 e- 20l / 3 (±e 9i -X) 

does not lie within an e-neighbourhood of any of the Stokes' lines, and z(h, A, —1), 
z(h,X,l) lie in different principal domains, one can always obtain (|i"4|) and (p!5|), 



and the asymptotics (|T3| ) will be valid. However, for A lying in the sector having 
its apex at T, and bounded by the rays 



T + re 29i/3 and 



T + r e 29i ' 3+2m / 3 for 6 < 
T + r e 2ei / 3 - 27ri /3 for 6» > 



r G [0, oo), it is easy to check that e~ 2km ^ 3 z(h, A, ±1) both lie in the same principal 
domain, for each k £ {—1, 0, 1}. Then it is also straightforward to check that as x 
ranges from —1 to 1, the function 

x^Re{(e- 2km/3 z(e,X,x)) 3/2 

which has been at the heart of our analysis, has a single maximum/minimum. 
Together with the identities (0), (f|) and (||), this means that there will be values 
of A such that equality holds in both ([L4]) and (|l|) - no matter what choices of j 
and k are made. Thus, (and without loss of generality assuming j = k = 0,) the 
set of A satisfying 

Re {(e-^V 1 " A)) 3/2 } = Re {(e" m / 3 (-e ei - A)) 3/2 } (16) 

lies in Y(V(—1),V(1)). We now examine this set in more detail. Expanding the 
Taylor series, we have 

/ e -2*/3, * _ A)) 3/2 = _ iA 3/2 e -« + 3 . A l/ 2 _ 3 A _ 1/2 « _ 1^-3/2^ + 

v v ;; 2 8 16 v ' 

and 

( e -2Wi/3(_ e « _ A)) 3/2 = _ zA 3/2 -« _ 3 l /2 _ 3 1/2 w i iA _3/2 2« + 0^-5/2) 

v v ;; 2 8 16 V ; 



as |A| — > oo. Dividing through by — i, this means that ( |16D will hold if and only if 

= 0(A- 5 ' 2 ) 



f8 

I 2 16 



as |A| — > oo. Putting A := pe^, this is equivalent to the requirement 

-(!)-£-(»-?)-<*^ 

as p — > oo. But then 

Im (A) = psin(0) 

= 2psinQ cos (fj 

= 2pcos(f|(-^sin (2e-¥)+0(p-*)\ 



2 1 I 24p 2 



o{ P - 1 ) 
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as p — > oo. By our definition of T, z(h, T, ±1) lies at the intersection of two Stokes' 
lines, and so 

Re{(^(/i,r,±l)) 3/2 } = 

showing that T certainly lies in the set of A satisfying (|T6|). Therefore, we deduce 
that Y(V(— 1), V^(l)) contains a curve from T asymptotic to the positive real-axis. 

Finally, we must examine what happens when z does lie on one of the Stokes' lines. 

Firstly, suppose Arg (z(h, A, 1)) = tc/3, corresponding to A lying on the ray centred 
at e dl and passing through V. Then k = 0, 1 are allowable, and one checks that if 
A lies on the segment [e 9 \ T), we have z(h, A, —1) G S_i. It follows by (H) that 

Re {(e- 2m/3 z(h, A, -1)) 3/2 } = Re {(z(h, A, -1)) 3/2 }, 

and so (|I1 ) and ( ^5|) cannot hold. However, if A lies on that part of the ray which 
extends past T (but not A = T itself,) then z(h, A, —1) G Si, and 

Re {(e- 2ni/3 z(h, A, -1)) 3/2 } = -Re {(z(h, A, -1)) 3/2 } 

causing (|14D and (|l^) to hold for j = 0, k = 1. 

An entirely similar argument holds when Arg (z(h, A, —1)) = vr, corresponding to 
A lying on the ray centred at — e dl and passing through V using (^), with j = 1 and 
fc = -l. 

Finally, the case where Arg (z(h, A, ±1)) = —n/3 is taken care of using (^), which 
shows that we may use allowable values —1 and to obtain fllil) and (|15). 

This completes the proof. 

In the case 9 = 7r/2; T = l/y/3 lies on the real-axis, and the figure Y(—i,i) has 
three linear 'arms'. When 6 = 0, the symmetric case, Y(—l, 1) is the semi-infinite 
interval [—1, oo), as is well-known from the theory of self-adjoint operators. 

We now give our main result. 



Theorem 2 Let 

H h f{x) := -h 2 ^^- + V(x)f(x) 



■ d 2 /(x ) 
dx 2 



act on L 2 (— 1, 1) with Dirichlet boundary conditions, where h > is small, and 
V(x) is the complex valued n-times piecewise linear function 



V(x) 



m\x + l\ xo < x < xi 
m<ix + I2 xi < x < xi 

m n x + I 



with —1 = Xq < X\ < . . . < x n = 1 and the rrii, li, i = 1, . . . , n complex constants. 
We assume for each i that if rriiXi + U = m i+ iXi + l i+ i, then 7^ rm + i. Put 
8i := Arg (mi), and, using our earlier notation 

n 

T:=\jY(V( Xi ),V(x i+1 )). 

i=l 



Let e > and N e Z + fre given. Then 

Spec(H h ) n {2 : |z| < iV} C Nhd(T; e) 
/or all small enough h > 0. 

Proof Our proof involves an analysis of the behaviour of the characteristic-determinant, 
i.e. the left-hand side of (M, as h — > 0. We give a proof for the case n = 3; the gen- 
eral case follows by a similar argument. For n = 3, the characteristic-determinant 
is given by 



un(-l) wia(-l) 

ttn(zi) m 12 (xi) -1*21 (aii) -1*22(^1) 

u' n (xi) u' 12 (xi) —u' 21 (xi) —u' 22 (xi) 

u 2 \{x 2 ) u 22 {x 2 ) -u 31 (x 2 ) -w 32 (x 2 ) 

u' 21 (x 2 ) u' 22 (x 2 ) -u' 31 (x 2 ) -u' 32 (x 2 ) 

1x31(1) u 32 (l) 



(17) 



and we must prove that for certain values of A G C, this determinant does not 
vanish as h — > 0. Expanding ([T7|), one obtains 

[(lill(-l)wi 2 (xi) - Ux 2 (-l)Uii(x 1 ))(u2 2 (x 1 )U2l(x2) ~ U 21 (xi)u 22 (x 2 ))x 

x(u 31 (l)u' 32 (x 2 ) -u' 31 (x 2 )u 32 (l))}- 

-{(lill(-l)Ml2(Xi) - Mi2(-l)M n (xi))(M 22 (Xi)M 21 (x 2 ) - u' 21 (xi)u' 22 (x 2 )) X 

x(w 31 (l)w 32 (x 2 ) -m 31 (x 2 )m 32 (1))}+ 
+ {(mh(-1)m / 12 (xi) - Ui2(-l)w / n(a:i))(w23(a!i)w2i( a: 2) - "21 (^1)1*22 (^2) )x 
x (1*31(1)1*32(^2) - 1/31(^2)1*32(1))}- 

— — 1)1*12(^1) - Ml2(-l)Mn(^l))(M 22 (xi)M 2 i(x 2 ) - M2l(^l)M 22 (x 2 ))x 

x (^1(1)1/32(^2) - u' 31 (x 2 )u 32 (l))} (18) 
where so far, no asymptotics are involved. 

Now, let e > and iV e Z + be as given in the statement of the theorem. Taking 
any 

\E{z:\z\< N}\Nhd(T;e), 

we can use the results of Proposition [L] to show that ( jI7[) is non-zero in the limit as 
h — *■ 0. Indeed, by the proof of Proposition [I], we can ensure that the asymptotic 
estimates 

wi 2 (-l)wn(xi) = o(wn(-l)wi 2 (xi)), 
«2i (^1)1*22 («2) = 0(1622(^1)1*21(^2)) 
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and 

hold, as h —>■ 0. Then, using the standard asymptotic expansions of the Airy 
functions 0, which give 

Ai(z) = ^ exp (-^ 3/2 ) (1 + 0(z-^)) (19) 



and 



Ai!{z) = exp (-\z* 12 ) (1 + 0(^ 3 / 2 )) (20) 



2^ 

as |z| — > oo, valid for all z such that |Arg (z)| < it; we see that, if 



z(h, A, Xj) := h 2/3 m i 2/3 ((miXi + U) - A), 



then 



d^ / 2 



exp 



dx2v / 7r 

dz 1/2 ^ 1/4 / 2 



(-^ 3/2 )(i + o(,~ 3 / 2 )) 



- — L '^=<>xp j -^ 3/2 ) (l + 0(^ 3 / 2 )) 



as | z | — > oo. Comparing this last expression with (|H]), and using / ~ g to mean 
that 

/(*) 1 



- 1 as /i — > 0, 



we obtain 



-|— Az(z(/i, A, Xi)) ~ — /i ^(mjX, + li) — X) 1 ! 2 Ai{z(h, A, Xi)) as /i — > 0. 
dx 

Moreover, similar calculations show that 

-^-Ai(e ±2m/3 z(h, A, ^)) ~ hrHOrmXi + k) - \)) 1/2 Ai(z(h, A, xA) ash^O. 
ax 

Reverting to our notation of (|3|), we will write 

u' a (xi) ~ h^cnix^unixi) etc. (21) 

as h — > 0, where it is important to note that the Cjj(xj), % = 1, . . . , (n — 1), j = 1, 2 
are independent of h. Then, since the constant terms Cjj(xj) are negligible in 
magnitude compared with the exponential terms h — > 0, the relations 

(pT|) imply that we also have the estimates 

ni 2 (-l)n' n (xi) = o(Mu(-l)M' 12 (ari)), 
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u' 21 (xi)U 22 (x 2 ) = 0(u' 22 (xi)u 2 l{x 2 )) 

u' 31 {x 2 )u 32 {\) = o{u 32 {x 2 )u 3 i{\)) 
u 21 (xi)u' 22 (x 2 ) = o(«2 2 (a;i)t4i(a:2)) 

and 

u 21 (xi)u' 22 (x 2 ) = o(u 22 (xi)u' 2l (x 2 )) 

as h — > 0. Returning to fll8|), we first use the above estimates (since we may ignore 
the sub-dominant term in each round-bracketed expression), and then the relations 
fl2lD again, to obtain the asymptotic estimate on the first of the curly-bracketed 
terms: 

{(Wll(-l)Wl 2 (^l) - U 1 2(-l)Uii(x 1 ))(u / 22 (xi)u 2 l(x2) ~ u' 2l (xi)« 22 (:E 2 )) X 
X («31 (1)m' 3 2 M - «31 (X 2 )U 32 (l)) } 

~ Wii(-l)ui 2 (a:i)t4!( a: i) u 2i(a!2)w3i(l)w / 32( a; 2) 

~ u n (-l)ni 2 (a;i)£: _1/2 C22(a;i)n 2 2(x 1 )u2i(x2)M3i(l)£ _1/2 C32(a;2)M32(a;2) 
= /i _2 [c 2 2(a;i)c32(a;2)](Mii(-l)Mi2(a;i)M22(a;i)M2i(x2)M3i(l)M32(x2)) 

as /i — > 0. Similar estimates apply to each of the remaining three terms in ( |I8|) i.e. 
{(wn(-l)wi 2 (xi) - ui 2 (-l)uii(xi))(u' 22 (xi)u 2l (x 2 ) - w 21 (xi)w 22 (a; 2 ))x 

x(n 31 (l)n3 2 (x 2 ) - %i(a; 2 )M3 2 (l))} 

~ /i~ 2 [c 2 2(a;i)c2l(x2)](Mii(-l)Mi2(xi)M22(a;i)M2l(a;2)M3l(l)M32(a;2)), 
{{u n (-l)u' 12 {xi) - U 1 2{-l)u' 1 ^X 1 )){u22{x\)u' 2X (x 2 ) ~ U 2 1 (xi)u 22 (x 2 )) X 

x{u 31 (l)u 32 (x 2 ) - u 31 (x 2 )w 32 (l))} 
~ /i _2 [ci 2 (xi)c 2 i(a;2)](«ii(-l)«i2(a;i)M22(a;i)«2i(x2)^3i (1)^32(^2)) 

and 

{(Mh(-1)m' 12 (Xi) -U 1 2{-l)u' 11 {x 1 )){u22{xi)U2l{x2) ~ U 2 i(xi)u 22 (x 2 ))x 

x(u 3 i(l)u 32 (x 2 ) - U 31 (x 2 )u 32 (l))] 

~ /i _2 [ci 2 (xi)c32(a;2)](ui 1 (-l)ni2(a;i)M22(a;i)«2i(^2)%i(l)%2(^2)) 

as h — > 0. Collecting these estimates together, we see that the characteristic 
determinant (|17|) tends asymptotically towards 

^ _2 {(c22(xi)-ci2(xi))(c 3 2(x2)-C2i(a;2))}(nii(-l)ni2(a;i)n22(a;i)M2i(a;2)M3i(l)M32(a;2)) 
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as h — > 0. The Airy functions Ai(z) have countably many negative real zeros, 0; 
and so by our choice of A outside Nhd(T; e) together with the proof of Proposition 
|I], we are assured that none of the Airy functions Uij(x{) vanishes. Therefore, the 
determinant (|T7D does not vanish in the limit as h — > 0, provided the 'constant' 
terms 

C22O1) 7^ ci 2 (^i) and c 32 (x 2 ) ^ c 21 (x 2 ). (22) 

Our choice of A ensures that each of the individual constant terms Qj(xj) is non- 
zero. Moreover, reviewing the proof of Proposition [l] and the identities (0)-©, we 
see that the choices for j and k are not uniquely determined. Therefore, it is always 
possible to ensure that (|22]) holds, even when V is continuous at some or all of the 
XiS. For example, if it happens that V(xi~) = V(xi+), then we choose j and k 
so that the constants Ci 2 (xi) and c 22 (xi) take different signs (by the calculations 
immediately above (|2"T|)). Thus, we deduce that such A cannot be an eigenvalue. 

It now just requires the following compactness argument to complete the proof. 
Let B(z; e) denote the open ball centred at z, with radius e. Our argument so far 
shows that for any 

Xe{zeC:\z\<N} 

such that 

B(\;e)nT = 

we have 

for all < h < E\, where E\ is some positive constant dependent upon A. Let 

M := {z eC : \z\ < N, and dist(z, T) > 2e}, 
so that M is compact. Then for all A G M 

B(X; e) fl Nhd(T; e) = 

and so 

MC |J B(\;e). 

AGM 

But by compactness this means that there exists a finite sub-covering 

n 

MC \jB(\ r -e Xr ). 

r=l 

Taking E to be min(£' Al , . . . , E Xn ) > 0, we deduce that for all < h < E we have 

Spec(H h ) n M = 
and this is equivalent to the statement of the theorem. 
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Remark 3 An important but subtle point, to note is that the zeros of 

(U1 1 (-1)U 12 (X 1 )U22(X 1 )U21(X2)U 31 (1)U 3 2(X2)) (23) 



as a function of A, are not the same as the zeros of (|17f) . However, by a similar 
argument to that of Shkalikov Q (i.e. using ([19])), one can readily show that 
along each of the bounded arms of the Y-shaped figures making up T, the zeros 
(eigenvalues) do converge as h — > to form a dense set. Finding an asymptotic 
expression for the density of spectral points along the infinite lines (in the direction 
of the positive real-axis) appears to be a much more difficult problem; and we have 
no results yet in that direction. 

To illustrate Theorem |2], in figure 3 we show a Matlab plot of the discretised version 
of the problem where the potential 



V(x) 



2ix + % for — 1 < x < 
(i + l)x for < x < 1 



We now return to the operator H$ h defined in the first section. 
Corollary 4 Let Hg h be the non-self-adjoint operator defined by 

H 5 , h := -h*^- 2 + V s (x) 
acting on L 2 (— 1, 1) with Dirichlet boundary conditions, h > , and 

V s (x) :-- 



i(x — 5) for x < 
i(x + 5) for x > 



with 5 > 0. Define S C C to be the double Y-shaped figure given by the line 
segments 

[i5,l/2Vs + i(l + 25)/2] 
[i(l + <y),l/2V3 + i(l + 2(J)/2] 



together with 
and 

together with 



[l/2v / 3 + *(l + 25)/2,+oo), 

[— i5, 1 /2>/3 — i(l + 2<5) /2] 
[-i{l + 6),l/2V3-i(l + 2S)/2] 



[l/2V3-i(l + 25)/2, +00). 
Then, given any e > and N G Z + ; we have 

Spec(H s ,h) H {z : \z\ < n} C Nhd(S;e) 

for small enough h > (see figure 1). 
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By analyticity, however, for fixed h > we have 

limSpec(if 5jh ) = Spec(if ,/i)- 

o— >0 

Hence, lini/^o lim^o Spe^if^) is contained within an arbitrarily small neighbour- 
hood of the line segments 

[i,l/V3],[-i,l/V3] and [l/Vs, oo) 

(see figure 2). Thus, the operations of taking limits do not commute, in the sense 
that as sets 

lim lim Spec(Hsh) nm nm Spec(iT,yh). 

5 Simultaneous limits for i?"^ 

In response to questions which we have been asked, we give an analysis for the 
situation in which 8 and h of Corollary [| are not independent of each other. 

Proposition 5 Defining the operator H$ h as above, and putting 

8 := h 1/p 

we have 

lim Spec(H h i/ P h ) = lim Spec(H h) if < p < 1 

and 

limSpec(if ft i /P ft ) = lim lim Spec(^) if p > 1. 

Proof Referring to @ and expanding, we see that the characteristic determinant 
of H s ,h is given by 

{(« u (-l)« 12 (0) - u 12 (-l)u n (0))(u' 22 (0)u 21 (l) - u 21 (0)u 22 (l))}- 

- {( Mll (-1K 2 (0) - ui 2 (-l)ui 1 (0))(u 22 (0)u 21 (l) - u 21 (0)u 22 (l))}, (24) 
whereas the characteristic determinant of Hq^ is given by 

u 21 (l)u 12 (-l)-u 22 (l)u u (-l). (25) 

Now, putting 

Ul2 (0) := Ai k (h~ 2 / 3 e~^ 3 (-i8 - A)) and w 22 (0) := Ai k {h- 2 ' 3 e-™l 3 {i8 - A)) 
it is clear by analyticity, that 

MO) ~ MO) and un(0)~« 21 (0) (26) 
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as 5 — > 0. Moreover, the calculations preceding ([H]) show that 

<4(0) ~ -5-f(-A) 1/2 « a (0) and ^ 2 (0) ~ <T p (-A) 1/2 Mi2 (0) (27) 

as 5 — > 0, for i = 1, 2. Therefore, using first (JZTJ) and then (|26|), the characteristic 
determinant fl24]) tends asymptotically towards 



2r P (-A) 1 / 2 (« 11 (-l)ttl2(0)«2l(0)M22(l)-Ml2(-l)«ll(0)«22(0)M 21 (l) 

~ 2r p (-A) 1 / 2 ( M22 (l)n 11 (-l) -%(l)«i 2 (-l)) 

as 5 — > 0. Then the zeros of ( p^D tend asymptotically toward the zeros of 
Rouche's theorem, explaining the behaviour of 

lim lim Sr>ec(Hs h). 



by 



Substituting 5 = h 1 '?, the character of lim^o Spec(if h i/ P h ) therefore depends upon 
the range of p for which 

uia(0) uu(0) 
as /i — > 0. Now, without loss of generality, and using our earlier notation, let 

^ 22 (0) . _ ^-i(^i) 
wi 2 (0) "~ Ai_ x {z 2 ) 

where 

Zx ■= h- 2/3 e m/3 (ih 1/p - A) and z 2 := h~ 2/3 e m/3 (-ih 1/p - A) 
so that, using the standard asymptotics (|6[) 

m 22 (0) 



Ml 2 (0) 



-1/4 / 

^ exp I 


:-u /2 ) 


(i + o(zT 


3/ 2)) 


-1/4 / 

2 2 exp I 






3/ 2)) 



£i 

2 2 

exp 



-1/4 

2 
~3 



exp 



3/2 3/2 
z l — ^ 2 



(i + o(^r 3/2 )) 



3/2 3/2 

1 ^9 



as h — > 0. But 



3/2 



3/2 

z 2 



hrW {{lb 1 /? - A) 3 / 2 - H/i 1 ^ - A) 3/2 | 

r l e ^/3(_ A )3/2 {(J _ ihVv/Xf/2 _ (1 + ^l/p/ A )3/2j 

/l -i e W3(_ A )3/2 1^ _ 3i/i i/P/2A + ...)" (1 + 3^ 1/p /2A 

/ l - 1 e 7ri / 3 (-A) 3 / 2 (-3z/i 1 / p /A + . . .) 
► 
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as h — > if and only if < p < 1. So, provided < p < 1 

«i 2 (0) 

and a similar calculation shows that we then also have 

«2l(0) . , 



1 as /i — > 0, 
«ii(0) 



as required. 
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